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We study superconductivity in the hadron-quark mixed phase of planar quantum chromodynamics 
(QCD) within the large N limit of a Gross-Neveu model modified by a repulsive vector term. At high 
densities, the combination of scalar attraction and repulsive space-like part of the vector interaction 
squeezes quarks into baryonic composite states, i.e., Dirac fermions with even numbers of bosonic 
vortices attached. The time-like vector component induces Cooper pairing between these Fermi 
surface modes. Remarkably, at zero temperature, competition between the quark density and mass 
destroys superconductivity via a Berezinskii-Kosterlitz-Thouless (BKT) phase transition driven by 
diverging chiral quantum fluctuations near criticality. Dissolution of logarithmically bound singlet 
diquarks is catalyzed by in-plane chiral mixing associated with Z 2 CZ 2 Z 2 chiral symmetry 
breaking of the Fermi surface into a transverse spin-polarized triplet ground state. We calculate 
the QCD phase diagram for quark chemical potential above the baryon mass based purely on Fermi 
surface considerations and find good agreement with results obtained by other methods. We address 
similarities between our quantum BKT transition and those found using holographic techniques. 


I. INTRODUCTION 

Presently, a thorough knowledge of the phases of quan¬ 
tum chromodynamics (QCD) is lacking and constitutes a 
significant hurdle towards our understanding of the rich 
phenomenology encoded within the standard model of 
particle physics m- Although QCD is notoriously dif¬ 
ficult to solve directly, much insight has been gained 
through a combination of laboratory experimentation, 
observational cosmology, and effective theoretical mod¬ 
els. For instance, data from heavy ion collision exper¬ 
iments tells us that for temperatures above the QCD 
scale, i.e., Aqcd ^ 200 MeV, and at low quark densi¬ 
ties the physics is described by a quark-gluon plasma. In 
contrast, we know considerably less about the physics at 
low temperatures and increasingly higher densities begin¬ 
ning near the baryon density, Aqqj^ ^ 1 fm“^. Although, 
it is true that at extremely high densities perturbative 
asymptotic freedom applies, as hadrons ultimately dis¬ 
solve into degenerate quark matter, asymptotic freedom 
does not apply for the intermediate regimes. Adding to 
the problem, first principle lattice simulations are im¬ 
practical here due to the infamous sign problem and 
empirical data is constrained by the fact that the re¬ 
quired quark densities are unreachable in the laboratory, 
occurring only in the deep interior of dense stelar ob¬ 
jects such as neutron stars. Nevertheless, various phe¬ 
nomenological methods such as the Nambu-Jona-Lasinio 
(NjL) gug and random matrix models mn\, among 
others, have lead to conjectured exotic phases including 
color superconductivity (CSC), color-flavor locking, and 
mixed phases characterized by the simultaneous coexis¬ 
tence of hadrons and quarks [ 2 ] . 

In this article, we investigate the interplay between 
chiral symmetry breaking (CSB) and color superconduc¬ 
tivity within the hadron-quark mixed phase of ( 2 + 1 )- 
dimensional QCD. We model the intermediate-to-large 
quark density regime of the QCD phase diagram within 


the large N semiclassical limit of the Gross-Neveu- 
Thirring model (see for example [ 8 ] and references 
therein). This model couples N species of fermions 
through scalar (Gross-Neveu) and vector (Thirring) in¬ 
teractions. As we are concerned primarily with interme¬ 
diate quark densities, a desirable feature of our model is 
the presence of terms that are equivalent to semiclassi¬ 
cal corrections coming from an external gauge field. This 
accounts for the early onset of asymptotic freedom associ¬ 
ated with perturbative QCD that applies towards higher 
densities. 

The character of the interactions in our model are as 
follows. The scalar term is taken to be attractive in order 
to allow for chiral condensation (qq) 7 ^ 0. In contrast, we 
have chosen a repulsive vector term. This is justified from 
muitipie iines of reasoning IMH- For exampie, both the 
instanton-anti-instanton moiecuie modei and the renor- 
maiization group approach resuit in repuisive vector cor¬ 
rections in their respective effective actions. From an 
empiricai standpoint, the equation of state for neutron 
stars requires a degree of stiffness that seems to be ac¬ 
counted for oniy by adding a repuisive vector term to 
standard NJL-type theories. The presence of a repuisive 
vector interaction piays a significant roie in the forma¬ 
tion of a coior superconducting phase. In particuiar, a 
space-iike metric signature produces a negative vector 
density interaction — (' 07 ^'^)^, which gives an attractive 
quark-quark channei and hence the possibiiity of diquark 
condensation (qq) ^ 0 associated with CSC. 

We wiii demonstrate how the interpiay between attrac¬ 
tion and repuision in our modei ieads to formation of the 
Fermi surface at iow temperatures by way of a Landau- 
Ginzburg-Wiison mechanism invoiving the ampiitude of 
the quark fieid. In essence, this is due to the repuisive 
space-iike part of the vector interaction ('07*'0)^. It is 
important to emphasize that in the fuii picture this oc¬ 
curs independent of, yet simuitaneous with, the forma¬ 
tion of biiinear condensates and hence chirai symmetry 
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breaking. Pauli blocking combined with the ensuing fi¬ 
nite expectation value for the quark amplitude forces the 
internal spinor phase to cycle around the circle in or¬ 
der to maintain a zero overall expectation for the quark 
field. Thus, we find that excitations of the Fermi sur¬ 
face are composite quarks: fermions with even numbers 
of bosonic vortices attached. The extra bosonic degrees 
of freedom effectively mitigate the repulsive exchange in¬ 
teraction between quarks. The net result is a large quark 
density stiffness consistent with what we know about the 
properties of dense QCD from neutron star considera¬ 
tions [ 18 ]. 

Our model exhibits a breakdown in the supercon¬ 
ducting phase through an infinite-order Berezinskii- 
Kosterlitz-Thouless (BKT) type phase transition. Sig¬ 
nificantly, this occurs as a conventional finite tempera¬ 
ture BKT transition along a boundary that intersects the 
chemical potential axis at zero temperature, hence sur¬ 
viving as a quantum BKT transition (QBKT). Detailed 
analysis reveals the QBKT mechanism to be fundamen¬ 
tally rooted in competition between the scalar part of the 
interaction and the time-like vector density. The former 
generates a quark mass which mimics an external mag¬ 
netic field perpendicular to the plane by introducing spin 
polarization asymmetry, while the latter generates a fi¬ 
nite quark density. The vortex degrees of freedom in our 
BKT transition are precisely the composite quarks that 
result from squeezing the quark amplitude and constitute 
the Cooper pairs near the Fermi surface. As we reduce 
the quark chemical potential (equivalently quark density) 
through the BKT transition the system favors quark- 
anti-quark over diquark condensation. Consequently, we 
find that diquark dissociation at BKT criticality occurs 
precisely when Cooper pairs transform from singlet to 
triplet states. 

The quantum BKT transition in (2 + l)d previously 
found for a D 3 /D 5 brane system [T 9 | using hologra¬ 
phy [20] [21] is of particular relevance to the results pre¬ 
sented in the present article. The initial holographic re¬ 
sults were subsequently extended to include various other 
driving parameters and systems [ 22 H 27 | . The QBKT 
transition in [ 19 ] was found to be driven by competition 
between applied magnetic field strength and quark den¬ 
sity, in contrast to standard thermally driven BKT. This 
scenario bears a remarkable resemblance to our problem 
wherein one finds QBKT to be driven by competition 
between the quark mass and chemical potential, where 
the effect of the former on the quark spin polarization is 
analogous to that of applied magnetic field and the latter 
directly proportional to the quark density. This parallel 
attests to the richness of four-fermion theories. Indeed, 
in this article we show that the four-fermion approach 
has the advantage of providing deeper insight by eluci¬ 
dating the microphysics underlying the BKT mechanism, 
in opposition to holography which defines the boundary 
field theory purely in terms of its bulk gravity dual at 
the unfortunate expense of a clear physical description 
for the BKT mechanism. 


This article is organized as follows. Section [IT] estab¬ 
lishes the foundation of our model by specifying the La- 
grangian density for attractive and repulsive four-fermion 
interactions. In Sec. |III[ we derive the structure of 
the Fermi surface based on low-temperature symmetry 
breaking, treating the quark mass as an external tunable 
parameter. In Sec. HYi we interpret excitations of the 
Fermi surface in terms of composite quarks. Here we un¬ 
cover the properties of quasi-quarks which result directly 
from the reduced compressibility of the system induced 
by quark-quark repulsion. In Sec. |V| we show that the 
dissociation of Cooper pairs in the superfluid state pro¬ 
ceeds via the Berezinskii-Kosterlitz-Thouless mechanism 


at both finite as well as zero temperature. In Sec. VI 


we 


connect to well-known established results by construct¬ 
ing, based on our model, the temperature-chemical po¬ 
tential phase diagram for QCD. In Sec. |vnl we conclude. 


II. THE MODEL 

We begin our analysis working directly from the micro¬ 
scopic dynamics of quark degrees of freedom interacting 
locally through Lorentz scalar-scalar and vector-vector 
terms. We compute the properties of the CSC diquark 
condensate from the resulting bound states of quarks of 
the same chirality. The Lagrangian density for our model 
is given by 


C = Co + Cs^Cv, (I) 


with the kinetic, scalar, and vector contributions given 
explicitly by 

N 

>Co = y] -m + /i7°) , (2) 

n=l 

= Y ’ ( 3 ) 


which includes explicit quark mass m and chemical po¬ 
tential jii. The superscript index refers to N species 
of quarks with scalar and vector couplings scaling like 
9 s^ 9 v ^ I/A^- This V-dependent scaling of the cou¬ 
plings will be important when we consider the semiclassi- 
cal regime. It follows that in the special case /i = m = 0 , 
Eq. 0 is symmetric under an SIJ{N)r C) SIJ{N)l chiral 
transformation. In the following we take the couplings to 
satisfy = Sgy = as there is no a priori relationship 
between these and this form leads to a spin-symmetric 
total interaction. Moreover, we will suppress the sum¬ 
mation over n, neglecting the species index, in which 
case Eq. 0 has an explicit discrete Z4 chiral symme¬ 
try. The covariant derivative — leA^ in Eq. 0 

couples the quark field to a non-dynamical gauge field. 
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as no associated field strength tr(F^jy)^ appears in the 
action. The reason for including this gauge potential will 
be made clear shortly. We work throughout in the two- 
dimensional Weyl formalism ^ '^2)^ using the 

first two Pauli matrices for the space-like part of our the¬ 
ory and the third one for the time-like part. 

An interesting feature of the model Eqs. 0-0 is that 
the large N limit of the time-like part of the vector in¬ 
teraction mimics weak coupling to a semiclassical gauge 
potential. To see this, consider a weak gauge interaction 
treated perturbatively and semi-classically by retaining 
only the scalar part of the four-vector potential. Ex¬ 
panded in a zeroth-order local static contribution plus a 
first-order nonlocal part, this four-vector potential reads 

^^(|r-r'|) [^o(|r-r'|), 0,0]^ , (5) 


with 


^o(lr-r'l) 


( 6 ) 


{l/a)Af^ (5(|r - r'|) + /(|r - r'|), 


where a is the fine structure constant. The first term pro¬ 
duces a point-like (mean-field) correction to the chemical 
potential and the second term a nonlocal correction. An 
effective chemical potential fl can then be defined as 


A(|r-r'|)^ 

fl + { 1 /q)A^q^ 


( 7 ) 

(5(|r - r'l) + {l/afA^Q'^ /(|r - r'|). 


But, the zeroth-order correction is identical to a Hartree- 
Eock term generated by condensation in the time-like vec¬ 
tor density Hence, coupling to a gauge field is ac¬ 

tually generated internally by the vector interaction and 
remains non-dynamical since the required field strength 
does not appear in C. Thus, promoting the kinetic term 
in Eq. Q to the status of covariance is warranted. 

The Lagrangian Eq. Q possesses a convenient symme¬ 
try, advantageous for modeling the hadron-quark mixed 
phase. Eirst, consider that one may view the parame¬ 
ter m as a dynamical scalar mass generated by a finite 
chiral condensate m = + mo, where mo is the 

bare mass. In this picture the chiral condensate is 

a Hartree-Eock mean field term generated by the attrac¬ 
tive scalar part of Eq. 0 . In a similar fashion, as ex¬ 
plained above, one may view the gauge potential Aq as a 
diquark condensate generated by the attractive time-like 
vector density through Aq = — 

with Aq^^ the bare potential associated with the diquark 
bound state. Eocusing on the density dependent parts, 
transforming between the quark-antiquark and diquark 
representations of >Cint, the terms —m and (g -|- Ao)j^ 
transform like 

(quark —antiquark) (diquark) 

-m; — 5 - , (8) 

(fi + Ao)y^: > -(V’V)- (9) 


Thus, the roles of the mass and chemical plus gauge 
potential are switched but their combination is invari¬ 
ant under transformations between the meson and quark 
matter representations. We will see that it is the com¬ 
petition between these two terms that is fundamentally 
responsible for dissolution of the diquark state. All of 
this strengthens the case for investigating Eqs. 0-0 as 
a plausible model for the coexistence phase of QCD near 
the superconducting transition. 


III. CHIRAL SYMMETRY BREAKING: 

FORMATION OF THE FERMI SURFACE 

The formation of a sharp Fermi surface in our model is 
closely related to competition between short range quark 
repulsive and attractive terms in Eq. Q. In this section 
we track the development of the Fermi surface at low 
temperatures by treating the quark mass as an external 
tunable parameter, revealing in the process chiral sym¬ 
metry breaking of the discrete Z4 symmetry of Eq. 0 . 
It is important to emphasize that in our treatment, chi¬ 
ral symmetry breaking is essentially generated within a 
Landau-Ginzburg-Wilson type framework for the quark 
field. We will see that at high densities the resulting fi¬ 
nite expectation for the quark amplitude forces stringent 
constraints on the internal phase of the quark spinor due 
to Pauli blocking for Fermi surface states. 

Formation of the Fermi surface is demonstrated by first 
expanding the spinor field in C in terms of single-particle 
Dirac states as 

[ci,i(t)xi,i(r), C2,i(t)x2,i(r)]^ > ( 10 ) 

and 

= Y AMxliir), C2,i{t)xli{r)] , (11) 

i 

where the summation index labels the single-particle 
states. We let the spinor functions ' 0 i(r) = 

[Xi,i(i*)5 X2,i(i*)]^ solve the time-independent massless 
Dirac equation 


{ia ■ V - eQ'tpiir) = 0 , (12) 

with eigenvalues e^. Using orthonormality f drxpi{r)xpj{r) 
= 6 ij, substituting the field expansions into the action, 
and performing the spatial integral, we obtain the action 
S = So + 5 'int with 


nhi 3 

So= dtYi^i + C2 i{t)dtC2,i{t) 

Jo i 

+ + €i\C 2 ,i{t)f 

+ [m- MT)] |ci,i(t )|2 - [m + A(T)] , 

/•ft /3 0^. 

5 'int= / dtY [ci,i(i)cxi(i)ci,fc(i)ci .;(0 

*^0 A A h. 1 






(13) 
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where the matrix for the spatially renormalized interac¬ 
tion strength is 

Sfjki - 9 ^ j ^’'X^,i(r)Xi,i(r)Xi,fc(r)xi,z(r), ( 14 ) 

with a similar definition for the lower spin component. 
Note that the temperature appears in the integration lim¬ 
its through P = {kBT)~^. With exponential forms for 
the coefficients Ci(2),i(^) = Ci(2),i exp(—at criti¬ 
cal temperatures Tc^i the combined quadratic and quar- 
tic terms develop local minima away from ci^i,C2,i = 
0, whose distances from the origin in spin space are 
eignevalue dependent. A significant point here is that, 
since we are interested in a finite quark density at low 
temperatures, the active modes are those lying nearest 
to the Fermi surface with energy ^ cf- The expansions 
in Eq. can then be approximated by sums over the 
single index . 

We should not be alarmed by the development of a 
finite expectation for the coefficients Ci(2), as these de¬ 
fine the overall and relative spinor amplitudes = 

\/c?|xiP +ci|x 2 p and -ci|x 2 p|, re- 

spectively, the latter encoding the spin polarization. 
Since Ci(2) are real, the phase structure of the fermion 
field p) resides in the space and time dependence of the 
decomposition Eqs. (p^-( 11 ), so that we correctly obtain 
= 0 even thoughTci^ 7^ 0 . We will this point in 
detail in Sec. hyi At the Eermi surface the Hessian ma¬ 
trix for Eq. ( [T^ can then be computed which yields four 
regimes defined by the following inequalities: 


1. for m + p — ej?< Q and —m + /i —ej ?>Q, two minima 

occur at C 2 =±v ^{ep — ft ci=0; 

2. for m + p — ei?>0 and —m-h/i — ej?<Q, tw o minima 
occur at C 2 = 0 , ci=±^/{{ep — ft 

3. for m p — €f<0 and —m p — ep < 0 , 
four minima occur at C 2 =^y^{p — m — ep)/g‘^^ 
Ci=±y {fi + m- €F)lg'^; 

4. for m + /i — ei?>0 and — m + p — ej?>Q , 

four minima occur at C 2 =^y^{fl m — ep) jg‘^^ 

Ci=±y(/i - m - eF)/g^. 


To delineate the temperature dependent regimes deter¬ 
mined by the above inequalities, we expand the chemical 
potential at low temperatures keeping only the first two 
terms. Eor a Eermi system /i(T) po—^2^^, where po is 
the chemical potential at T =0 and the coefficient /i2 can 
be determined by standard methods. The broken sym¬ 
metry phase corresponding to four minima in Eq. ( 13 ) 

is delineated by the condition T < Tc^\m)^ where the 
mass-dependent critical temperature is 


T^im) = y(l/M2) [mo - (to + ef)] , ( 15 ) 

and we have defined po = — A^. In contrast, the 

phase with two minima is determined by the condition 
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FIG. 1: Chiral symmetry breaking assoeiated with formation 
of the Fermi surfaee. (a) The temperature-mass plane phases 
in the development of the Fermi surface are separated by 
Landau-Ginzburg-Wilson type transitions. Ghiral symmetry 
of the full theory in each region is indicated as either broken 
(x^^) or retained (x^). In regions I and V (T > no well- 

defined Fermi surface exists. In region II (Ti^^ < T < 
the Fermi surface is formed in the presence of a large quark 
mass or a small mass at moderately low temperatures. Re¬ 
gion III (T < Tc^\ m 7 ^ 0) corresponds to a well-defined 
Fermi surface in the presence of a small but non-zero quark 
mass. In region IV (T < Tc^\ m = 0), the Fermi sea is 
comprised of massless quarks, (b) Minima (ground states) of 
the effective potential in Eq. ( p^ are depicted as disks in the 
spin-space coordinate plane for each region in (a) with spin 
up and down along the vertical and horizontal axes, respec¬ 
tively. The inherited chiral symmetry (degeneracy) group of 
the Fermi surface (FS) is also indicated. 


Tc^\m) <T < Tc^\m)^ where 


(1), 


rib(m) = 7(l/M2)[/io + (TO-eF)]. 


( 16 ) 


Figure depicts the various temperature-dependent 
regimes involved in forming a sharp Fermi surface. Note 
in particular the quantum critical point rric = \eF — Ao|- 
A key point in our analysis is that the critical tem¬ 
peratures, Eqs. ( 15 )-(H^, depend on the Eermi energy 


ep- This is because we have only considered modes for 
the expansions in Eq. ( 13 ) with energies near ep- Higher 


temperatures blur the Eermi surface allowing for active 
modes deep within the Eermi sea. When this occurs. 
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the structure underlying Eqs. (15)-(16) is smoothed out 
in that we can no longer identify a single effective po¬ 
tential for the action Eq. (13). Instead, what we see at 
higher temperatures is a quantum blurred band of “criti¬ 
cal” temperatures with no particular energy (or temper¬ 
ature) single-out as special. A sharp transition in the 
Eermi sea occurs only when a particular eignevalue dom¬ 
inates the expansions in Eq. (13). 




IV. EXCITATIONS OF THE FERMI SURFACE 

The development of finite expectation for the quark 
amplitude requires a completely random internal phase in 
order to maintain the required fermion field condition (t/^) 
= 0. Consider a general two-spinor field parametrized as 
'ip = exp{i0) [exp(i0/2), exp(—i0/2)] , where for single 
particle states the external phase angle 0 = Op = p ■ r/h 
is related to the momentum of the particle and the inter¬ 
nal geometric phase p = p = tan“^(p^/pa,) is identified 
with the polar angle and determines the direction of mo¬ 
mentum flow in the plane. Note that p provides the cor¬ 
rect half-angle property characteristic of a fermion field. 
Hence, p parametrizes a degenerate manifold for a given 
value of |p|. In order to satisfy the condition {pj) = 0, 
for a given energy eigenstate pj must form a superposi¬ 
tion over all possible momentum directions, i.e., values 
of p must be randomly distributed wrapping around the 
polar angle p. The simplest way to satisfy this is for p to 
he within the range 0 < p < 27r(2n + 1), where n G Z+. 
The contribution from n in the upper bound describes 
2n bosonic vortices attached to the quark field with the 
Fermion geometric phase contained in the second term. 
Thus, the Fermi surface is in essence described by two 
orthogonal sets of momentum direction distributions. At 
each point in space these can be chosen to be the radi¬ 
ally diverging (hedgehog) and rotational (vortex) flows, 
as depicted in Fig. |^a). In each case the momentum 
angle wraps around the polar angle 2n times. 


(b) 



FIG. 2: (color online) Quasiparticle structure of the Gross- 
Neveu model with repulsive vector interaction, (a) Local basis 
for single particle states comprised of superpositions of radial 
and vortical momentum states. The squeezed quark ampli¬ 
tude induces superpositions of different momentum directions 
at each point associated with a random internal phase of the 
spinor field. This effective introduction of topological order is 
required in order to maintain the fermion condition {fj) — 0. 
The extra phase winding describes an even number of bosonic 
vortices (orbital angular momentum L) attached to a single 
quark (spin S) creating a composite fermion, reminiscent of 
quantum Hall systems, (b) Density fluctuations over a uni¬ 
form background expose a net extra phase winding, as adja¬ 
cent fluid currents no longer cancel. Bound regions of surplus 
(light) or deficiency (dark) of fluid appear as oppositely ro¬ 
tating vortices. Vortex unbinding takes place when m = ttIc, 
at which point S becomes polarized transverse to the plane, 
in the direction of L. Note that the radial basis in (a) is 
associated with non-equilibrium time-dependent processes. 


A. Composite Quarks 

Let us now see how quasiparticles emerge in this pic¬ 
ture. Consider a density fluctuation over a uniform back¬ 
ground in the vortical basis at the Fermi surface, depicted 
in Fig. [^b). The resulting density gradient leaves some 
of the background vorticity “exposed” since contributions 
from neighboring fluid elements no longer cancel. As a 
result, a net current flows in the regions with elevated 
or depressed density, circulating in opposite directions in 
each case. This localized circulation reflects the statis¬ 
tics of the surrounding medium. That is, quasiparticles 
are fermions with bosonic vortices attached. This can 
be understood by recalling that the extra attached vor¬ 
tices were required in order to maintain overall Fermi 
statistics as the quark amplitude develops a finite ex¬ 
pectation value, effectively forcing quarks into a densely 


packed squeezed state. Under such conditions the at¬ 
tached vortices provide bosonic screening between nearby 
fermions. The significance of this point can be better 
appreciated by recalling that fermions and bosons ex¬ 
perience, respectively, repulsive and attractive exchange 
interactions. Hence, high density quarks are sparred the 
full impact of exchange repulsion by forming composite 
structures comprised of a bosonic vortex with a fermionic 
core, a more energetically favorable configuration than 
the bare fermion. 


[1] The phenomenon of composite particles is quite familiar in quan¬ 
tum Hall systems. For example, see Ref. m for a thorough 
treatment of composite fermions, and Ref. m for an interesting 
overview of general composite particles in quantum droplets. 
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B. Characterization of Quasi-quarks 


In order to gain a more quantitative understanding of 
Fermi surface collective excitations, we proceed next to 
work in the hydrodynamic formalism through the canoni¬ 
cal Madelung decomposition of the fermion field in terms 
of Weyl component amplitude and phase fields 




cos7(r, t), sin7(r, t) , 


(17) 


where 'd(r, t) and ip{T, t) are the overall and relative 
phases, respectively, with the local spin-up (down) densi¬ 
ties pi( 2 )(r, t) related to the quark amplitude ^(r, t) and 
chiral angle 7(r, t) by pi = and p 2 = r^^sin^y. 

We will soon see that this choice of spinor decomposi¬ 
tion is crucial in order to identify the underlying mech¬ 
anism of the QBKT transition. For convenience we now 
shift to using the total energy functional associated with 
our Lagrangian density. Incorporating the decomposition 
Eq. ( pT] ), then, after substantial algebraic manipulation, 
we extract the total energy for the Weyl spinor compo¬ 
nent fields 


E = dr 


h 


i sin27 • V'd + - sin27 • Vp^ (18) 


+ 7 ^cos 27 n^ • Vy + ip‘^ |n^ x Vyl 


-h p^ sin 27 • V(^ sin 27 |n^ x V(p\ 

+ 7^ n^M -1- ^ (cos^y + sin^y) 


Note that at very low temperatures, Eq. (18) is approx¬ 
imately equal to the Landau free energy (n^lecting the 
entropy term). Erom here forward we will work in this 
time-independent formalism where covariance is not ex¬ 
plicit but which favors a more physically intuitive per¬ 
spective. 


The gradient terms in Eq. (18) appear as scalar and 


vector products of the unit vector = (cos(/9, simp) 
with gradients of the various fields. These terms can 
be interpreted geometrically since the angle (p describes 
the direction of the Dirac current. Hence, is a 
unit vector along that direction, which identifies the 
gradient terms in Eq. (18) as directional flows for the 
fields p and y with respect to the quark momentum 
^ V'd. The mass and interaction terms in the last line 
of Eq. (18) depend on both the overall density p‘^ as well 


as the chiral field y through the spin polarization vector 
= (cosy, siny) and in the quartic interaction. The 
temperature dependent 2x2 mass matrix is defined as 
M(T) = diag [m — jl{T ), —m — /i(T)], where we have in¬ 
dicated the chemical potential to be explicitly dependent 
on the temperature T. 

To reveal the quantum BKT transition in our model, 
we focus on the broken chiral phase of the Eermi surface 
in Region III of Eig.[^(T ^ Ti^^) near the critical point 
m < rric and expand the massive fields near a minimum 


of the energy Eq. (18). This yields the energy functional 


E = Eq + Eint for the field fluctuations 5p^ (5y, 'd, and (p 
given by 


Ec 


— dr^s^uiSp^\\dp T 


(19) 


1 


+ (cmV||(5y + iV±6j) + 

2 ^rj^m (zV||<y?-V_L (/?)], 
-Eint = Jdr [8|toc| {Sp + Sm2\mc\5p/g^) (20) 

- 4:Sm\m - mc\Sp , 
where the coefficients are 




( 21 ) 


Om — 


— 


(1-2^7^) 


1^1 

\mr\ 


— 2^7 -^—r\m + m — rric 72, (22) 

I'^cl 

(1 — 2(5y^) -^—\m + — md^^^ 


-2^7 ■ 


rric 


(23) 


Expanding about a minimum of the energy shows the 
emergence of two massive fluctuations, Sp and 6^ in 
Eq. (20), associated with radial and rotational directions 
in spin space, and two massless modes d and p. Classi¬ 
fying the fluctuations in our system, we find (parallel || 
and perpendicular T notation referring to directions rela¬ 
tive to the quark current): 1) a longitudinal spin density 
wave, or density compression wave, 6p parallel to the gra¬ 
dient of the quark phase'd; 2) a massive spin-wave fluc¬ 
tuation in the local chirality encoded in the chiral fluctu¬ 
ation 3) a phase fluctuation, or spin wave, associated 
with the quark current V'd; and 4) fluctuations p in the 
direction of the quark current. Encoded in V(^ are the 
longitudinal “snake” mode \/\\p, which nucleates quan¬ 
tum turbulent flow in the quark bilinear condensates, and 
a transverse spin wave \/±p. This last mode introduces 
spatial divergence into the quark current, inducing local 
fluctuations in the quark density. 

The gradient terms in Eq. (19) may be cast into a more 
enlightening form using geometric arguments. Specifi¬ 
cally, fluctuations in the internal angle may be shown to 
convert to ■ Vp = \^\\p\ |V x n^| and |n^ x V^| = 
|V' 11 ,^ 1 ; the former encapsulates orbital and 
spin currents; the latter accounts for time-dependent con¬ 
vection flow (see m for the Gordon decomposition of 
Dirac current and [SD Eg for the Weyl formulation of 


[2] In the full dynamics associated with Eqs. ( |19| >-( [20} the parallel 
and traverse gradients in the phase cp are coupled and describe 
time-dependent fluctuations in the quark momentum held V-d 
akin to electromagnetic held propagation. 
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Maxwell’s equations). To relate the ip fluctuations to the 
divergence and curl of the direction field one expands 
the two orthogonal directional derivatives of p as follows 


V_L(^ = (7^7r/2nv) ■ (24) 

d d 

= —^\iip—p^co^p — p (25) 

ox oy 

= V • (cos(/ 9 , sinp) (26) 

= V-n^, (27) 


and 


V||^ 


• \/p 

d 


d 


cosp—p + smp — p 
ox Oy 

|V X (cos(/ 9 , sinp) I 

|V X n^l . 


(28) 

(29) 

(30) 

(31) 


Note that we have used the 2 D rotation matrix 71^/2 to 
rotate the direction field by 90 degrees. Invoking the 
substitution |V^(^| ^ |V • n^|, |V||(/?| ^ |V x n^| 

provides a more physically intuitive perspective. In¬ 
deed, as the field describes the direction of momen¬ 
tum flow, the divergence and curl forms here reveal pre¬ 
cisely t he sam e hedgehog and vortex structure presented 
in Sec. |IV A and illustrated in Fig. [^a). We emphasize 
again that in the hydrodynamic picture, the curl of the 
direction field embodies the fermion spin and orbital 
angular momentum as well as additional quantized rota¬ 
tion of attached vortices, whereas the divergence of 
describes the convection part of the fermion current into 
or out of a volume element. 


incorporate the fundamental quark spin (Pauli vortex) 
and the extra attached bosonic vortices, both contained 
in the spin/orbital term of the Gordon decomposition 
of the Dirac current. By casting our problem in terms 
of a gas of such defects and integrating out the massive 
density and chiral modes, 6r] and we will find that 
fluctuations in the former produce a hard-core delta func¬ 
tion repulsion, whereas fluctuations in the latter provide 
a pairwise logarithmic attraction between defects whose 
coefficient vanishes at the critical point m = rric. 

At ultra-low temperatures around a chosen ground 
state, to a good approximation, massive fluctuations can 
be treated as non-dynamical fields with each propagator 
contributing an overall factor of the inverse mass. Near 
a quantum phase transition, however, th e chiral mode 
5j becomes massless decreasing as ^ ^ymc{mc — m), 
in which case the corresponding propagator contributes 
a logarithmically divergent energy. Thus, at zero tem¬ 
perature the diquark dissociation energy at the critical 
point m = rUc is supplied by long-wavelength quantum 
fluctuations in the chiral field. 

The general structure described here is precisely that 
of a 2 D Coulomb gas with an additional infinite repul¬ 
sion at the defect cores. Consequently, implementing the 
standard renormalization group analysis reveals the BKT 
signature exponential decay of quasi-long-range order for 
melting (when T 7 ^ 0) or dissolving (for T = 0) of the 
diquark condensate. The crucial point to emphasize here 
is that, in order to work, the precise mathematics that 
maps composite quarks to the Coulomb gas degrees of 
freedom requires a relativistic Dirac structure, as the lat¬ 
ter is endowed with a peculiar chiral structure absent in 
ordinary BCS theory. 


V. ZERO-TEMPERATURE 
BEREZINSKII-KOSTERLITZ-THOULESS PHASE 
TRANSITION 


A. Derivation of the QBKT Transition 


Since its inception, the BKT mechanism has found 
application in a wide range of physical systems char¬ 
acterized by a transition from algebraic to exponen¬ 
tial behavior of long-range order j33ll35] . The topo¬ 
logical nature of the BKT transition circumvents the 
well known Mermin-Wagner-Hohenberg no-go theorem 
forbidding finite-temperature spontaneous breaking of a 
continuous symmetry in dimensions < 2 [36ll38] . Al¬ 
though in 2D spontaneous symmetry breaking at zero 
temperature is still possible. In the particular case of the 
holographic BKT transition for a D3/D5 brane system 
presented in m, the underlying mechanism is rooted 
in the violation of the Breitenlohner-Freedman bound on 
the gravity side and the fact that parameters such as den¬ 
sity and magnetic field have the same mass dimension. 

In this section we derive the contribution to correla¬ 
tions in the superconducting phase of QCD coming from 
defects in the quark field, i.e., the composite quarks in¬ 
troduced in Sec. Eva] As stated earlier, these defects 


The effective energy for bound defects is obtained by 
Gaussian integration over the massive modes Sy and ^7 in 
two steps. First by integrating over the quadratic contri¬ 
bution from the Sy modes in Eqs. (p^-(2Q), then gather¬ 


ing the ^7 quadratic terms in the resulting expression and 
integrating over these. At each step the functional inte¬ 
gral is performed by Fourier transforming to momentum 
space, then applying the standard Gaussian prescription 
(saddle-point approximation), re-exponentiating the re¬ 
sulting determinant, expanding about the extremum of 
the action, and finally re-expressing the result in real 
space. This procedure gives an effective theory for the 
spin wave (smooth) component of the quark current and 
a quark spin (defect) contribution from the vortical term. 
The defect contribution to the effective energy near crit¬ 
icality {m <mc) is 


E =Jdrdr'ip^{r')[pr,6{r - r') + p^\n (r - r')]ipd.{r), (32) 









with 


Pr] — 


64m^ 


(m + mc){mc — tu) , P^ = 


Airm. 


Pr] 5 ( 33 ) 


and where the defect field is given hy ipd = z • (V x 
n^^). The logarithmic contribution in the second term 
of Eq. (32) is key, and comes from the long distance be¬ 


havior in the momentum integration of the 5j Green’s 
function. Defect fields of the same chirality (opposite 
spin and momentum) are now coupled through a repul¬ 
sive contact term with coefficient which comes from 
integrating out the amplitude field, and an attractive log¬ 
arithmic term with coefficient coming from the chiral 
field. Note that both coefficients in Eq. (33) vanish at 
m = rric but at different rates. Diquark binding is de¬ 
picted in Eig. I^a). 

Contrasting this with the result far from criticality 
(m ^ rUc), we find that the second term in Eq. ( [^ is 
an attractive delta function in the separation ^(|r — r'|), 
which reduces the strength of the hard-core repulsion 
coming from the first term. This happens since the mass 
of 6y fluctuations in this regime is large compared to the 
characteristic low-temperature momentum. Thus, in ad¬ 
dition to the unbinding that occurs for m > rric^ a cross¬ 
over occurs between two regimes: logarithmic binding 
around m < rric to unbound defects for m rric. 


Equation (32) can be mapped to the 2 D Coulomb gas 


as follows. Note first that for a closed path encircling 
delta function defect sources we have 

^ ' ds = j (drz) • V x (34) 

^ V X n^(r) = 27rz^n^ (5^(r - r^) , (35) 


where ds is a tangent differential vector along the path. 
“Magnetic” charges for delta function defects at positions 
r^, Tj, are denoted as rP-. Thus, in terms of an analog 
magnetic held notation V x = V x = zV^M^, 
we then obtain V^M(i(r) = tt ^^(r — r^), which has 

the solution Md(r) = tt z In (r — r^). Incorporating 

these results into Eq. ( [^ and performing the integration 
over r and r' gives 

-E = TT^ y] <n}’ [Pn - r^) + p^\n{Yi - r^)] , (36) 

ip 


which leads to the diquark condensate correlation length 
obtained through standard renormalization techniques 


exp 


c/y\m - mc{T)\ 


, where c is a constant and 


we have included finite temperature effects. The correla¬ 
tion length diverges exponentially near the critical point 
mc(T), consistent with the BKT theory: at T = 0, a 
QBKT transition occurs at m = mc(0); for T 7 ^ 0, mc{T) 
gives the critical point for a standard BKT transition, 
modified by quantum mechanical corrections. 

The zero-temperature transition displayed in Eq. (36) 
is driven by the quantity p-\- Aq — m (see Eig. [^, which 


multiplies the logarithmic term and relates to the mass of 
chiral fluctuations Sy. The remarkable similarity of the 
second term in Eq. (36) to the free energy in conventional 


BKT for a system of size R comprised of vortices with 
radius a given by 


F = E-TS = {n- 2kBT) \n{R/a) 


(37) 


is evident by noting that the sum of chemical potential 
and vector density in our system measures the energy 
of an isolated defect, analogous to k in Eq. (37), with 


the scalar mass analogous to the entropy factor 2k bT. 
Hence, the “thermal” energy required to overcome vor¬ 
tex binding in the QBKT transition is evidently supplied 
by quantum fluctuations as m approaches criticality. In 
Sec. |V C[ we will see that these fluctuations correspond 
to induced admixing of triplet pairs into the pure singlet 
Cooper pair states. 


B. Analogy to 2D Magnetic Spin System 


As we have shown, at zero temperature the Eermi sur¬ 
face can be classified according to the symmetry of the 
corresponding large N semiclassical theory: formation 
of the Eermi surface breaks the chiral symmetry of the 
full finite temperature theory due to the sign change of 
the quadratic terms in Eq. (13). Consequently, the quark 
spin acquires a well-defined polarization transverse to the 
plane analogous to the quantum spin problem in the pres¬ 
ence of a transverse magnetic held at low temperatures. 

In fact, the analogy is made evident by relating our 
parameters to a fictitious magnetic held defined as S = 
(S||, S±)^, where B\\ = —{p + Aq) and Bb = —as 
shown in Eig. [^b). Under chiral symmetry breaking, 
the symmetry of the initial finite-temperature theory is 
inherited by the degeneracy of the ground state, i.e., the 
Eermi surface. In the case where the scalar mass m = 0 
(S_L = 0), the Eermi surface has a Z 4 degeneracy, related 
to the number of minima in the effective potential for the 
quark held in Eq. (13). This describes chiral symmetry 
of a Eermi surface comprised of two left-handed and two 
right-handed chiral ground states with the quark spin po¬ 
larized along the plane. Tuning the mass away from zero, 
but well below the critical value 0 < m < rUc (analogous 


to rotating B out of the plane), breaks the Eermi surface 
degeneracy to Z 2 C Z 2 . Here, what is initially a balanced 
superposition of transverse spin polarizations, equivalent 
to a well-defined in-plane chirality, is subsequently bro¬ 
ken as the quark held is driven towards an admixture of 
left and right chiral states. When m > rric = +Aq, a fi¬ 

nal round of chiral symmetry breaking occurs to a ground 
state with Z 2 degeneracy associated with spin polarized 
in the positive z direction, i.e., an equal superposition of 
left and right chirality for the quark held. 

Thus, at T = 0 tuning the ratio m/(/i+ Ao) (or B^/B\\ 
in our analogy) through unity, dissolves the diquark con¬ 
densate through a mechanism analogous to magnetic 
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(a) 




Pauli matrix as acting externally as a discrete symmetry 
transformation ^ 


The Madelung decomposition of the Weyl spinor field 
takes the form exp(i6>)[cosy, ±exp(i(/)) siny]^, 

where the various parameters are functions of the co¬ 
ordinates. Here, as discussed, the additional field y is 
generally massive and therefore non-dynamical at zero 
temperature (except near a phase transition). Its ex¬ 
pectation value is related to the particular ground state 
of the system through the average polarizat ion vector 
P = (cos(y), sin(y))^ oc (ci, 02)^ from Sec. IIP where 
proportionality accounts for a normalization constant. 
For illustration consider that we are in the = (1, 1)^ 
ground state of the manifold with Z4 symmetry, which 
corresponds to (y) = 7r/4 in the right hand chiral ground 
state. Consider then increasing the quark mass by a 
small amount so that the ground state symmetry is bro¬ 
ken Z4 ^ Z2 (8) Z2. To expand about the new minimum 
we take (y) ^ 7r/4 + (5y, so that 


FIG. 3: (color online) Diquark in the Gross-Neveu model with 
repulsive vector interaction, (a) The binding force between 
composite quarks arrises from the quark chemical potential 
jjL and a mean field Hartree contribution from the time-like 
component of the vector interaction ^o- The scalar mass 
m competes with these by adding repulsion between quarks, 
favoring chiral condensation instead. When m — uic = 
/X + To, unbinding of composite quarks is driven by a quan¬ 
tum Berezinskii-Kosterlitz-Thouless transition, (b) Binding 
of composite quarks cast in terms of the analog magnetic field 
B = (i3||, Bi_)^, where B\\ — —{fi + To) and B± = —m. For 
\B^\ < \B\\\ the quark spin polarization is mixed between 
in-plane and transverse polarization. The QBKT transition 
occurs for B^ = B\\. Note that the quark spin polarization is 
purely transverse for \Bi_\ > \B\\\. In this sense, one can see 
that unbinding of composite quarks is analogous to conven¬ 
tional magnetic catalysis in systems with an applied external 
magnetic field. 


catalysis, which has been studied extensively for rela¬ 
tivistic fermions in the presence of an external applied 
magnetic field (see for example Ref. [39] and references 
therein). 


C. Chiral Fluctuations 

The quantum phase transition at m = rric is driven by 
diverging quantum fluctuations in the chiral field Jy. To 
explain precisely what we mean by “chiral field” in the 
present context consider first that in (2+l)-dimensions 
the right and left massless chiral spinors can be expressed 
a-s i^R^L = p^/^exp(i^) [1, ±exp(i0)]^, with the positive 
(negative) sign associated with the right (left) handed 
spinor, and 0 = tan“^(p^/pa,). In 2D, the group SU(2) 
gets reduced to U(l) (8) Spin(2) with only two of the Pauli 
matrices <71^2 playing a dynamical role. The equations 
for left and right hand spinors decouple with the third 


P ^ [cos(7r/4 + (5y), sin(7r/4 + (5y)]^ (38) 

« (1, If+ 57 ( 1 ,-1)^ (39) 

= pf+<57Pf. (40) 


Hence a small fluctuation 6^ around a particular ground 
state, here introduces a component of the opposite 

chiral ground state P^^^ Stated in words, the fluctuation 
(5y describes quantum tunneling into an adjacent ground 
state as shown in Fig. 


Functional integration over the massive fields couples 
quarks of the same chirality, i.e., states of opposite mo¬ 
mentum and spin polarization, forming singlet Cooper 
pairs. Chiral fluctuations in the second term of Eq. (40) 
affect the overall spin of these Cooper pairs. Consider 
the diquark state { 2 p]^ 2 pR). The second term in Eq. (40) 
induces the fluctuations 


, ( 41 ) 

from which one sees that first-order corrections in the 
chiral field 6^ comprise spin-triplet states (unstable in 
ordinary BCS theory) whereas the second-order term de¬ 
scribes tunneling into the opposite singlet state. Thus, 
the dominant effect of the field is to drive the indi¬ 
vidual chiral quarks that make up the diquark towards 
mixed chiral states. By tuning m rUc, fluctuations 
in (5y become more pronounced as adjacent minima in 
the effective potential are brought closer together (see 
Fig-i, reducing the height of the barrier between them 
and mixing the opposite chiral state into the quark field. 
Einally, the quark field becomes an equal admixture of 
left and right chiral states at m = rUc- Erom this discus¬ 
sion, unbinding of composite quarks at criticality is then 
seen to result from the instability of the triplet bound 
state. 
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FIG. 4: (color online) Quantum tunneling between adjaeent 
ground states of the effeetive potential. Tunneling in spin 
space is depicted from the to the ground state driven 
by the S'y field. As the quark mass is tuned m rUc the two 
potential minima approach each other (see depiction of region 
III in Fig. Bb)) significantly increasing the tunneling rate. 
Consequently, near criticality, Cooper pairs undergo extreme 
fluctuations between the singlet and triplet bound state until 
finally dissociating at m = ttIc. 


VI. QCD PHASE DIAGRAM 

In the final part of our analysis we translate our re¬ 
sults into the more familiar form in terms of the quark 
chemical potential /i, which can be related directly to 
the baryon chemical potential or the quark density. The 
region of the QCD phase diagram that emerges is the 
transition region between hadronic matter, defined by 
(qq) > (gg), and quark matter, (gg) < (gg), which 
is sometimes referred to as the coexistence phase of 
QCD [2j [3]. This is intimately related to the idea of 
hadron-quark continuity, i.e., a smooth transition from 
superfluid/superconducting hadronic matter to super¬ 
conducting quark matter. 

Another way to view this region is in terms of what 
has been dubbed “quarkyonic” matter [40], wherein the 
fermionic sector is defined by a mixed system comprised 
of free quarks as well as baryons. Unique features of 
this phase include the possibility of chiral symmetry 
restoration in the presence of confinement. We show this 
to be true for our model: chiral restoration occurs at 
low temperatures and for large chemical potential. Our 
model also exhibits the effect of a baryonic “skin” at the 
Fermi surface, another characteristic of the quarkyonic 
phase characterized by a layer of baryons of thickness 
^ Aqcd measured from the Fermi surface that interpo¬ 
lates smoothly into the pure quark phase residing deep 
within the Fermi sea. 

As we have seen, this baryonic picture of the Fermi 
surface is realized in our model through the composite 
quasi-quarks derived in Sec. EVA) the lowest of which 
are naturally identified with the lightest baryon states. In 
extracting the QCD phase diagram from our model, we 
must keep in mind that a Fermi sea can only form when 
the chemical potential exceeds the baryon mass M. This 


reasoning is crucial for identifying a lower energy bound 
for the viability of our model. Hence, from our analysis 
of the Fermi surface we find that the tri-critical point for 
the onset of superconductivity occurs at /i = M. 

To construct the temperature-chemical potential phase 
diagram, we recast the boundary conditions from Sec. |m] 
for the various chiral broken phases of the Fermi sur¬ 
face using a quadratic approximation for the scalar mass 
as a function of temperature m(T) = mo — m2and 
taking the chemical potential to be the free parameter. 
We then re-express the effective chemical potential in 
terms of the gauge potential ft = ja — Aq, and include 
the dynamical contribution to the gauge field through 
^0 = + Aq^\ with > 0. This gives 


T A 


X 

deconfined phase 
(QGP) 


V mo/ m 2 Z 2 hadronic'‘•C jU 

phase 

[qq) ^ 0 : 

{qq) = 0 : QBKT 



{qq) ^0 / 

(qq) ^ 0 / ^4X' 

I restorati 


m > 1 


M -\-mo 


■0 


/i 


FIG. 5: Temperature-ehemieal potential phase diagram for 
the Gross-Neveu model with repulsive veetor term. C hira l 
symmetry breaking phases are indicated based on Eqs. (42)- 
(43) associated with formation of the Fermi surface. The Z 2 - 
symmetric hadronic phase with vanishing diquark condensate 
lies to the left of the QBKT transition line. Just to the right 
of the QBKT line lies the superconducting phase described by 
a Coulomb gas of Cooper pairs with Z 2 C Z 2 symmetry. For 
very large fi the system crosses over to free quarks indicated 
by the dashed curve to the lower right, obtained by shifting 
the lower curve in Fig. |III| to smaller values of m. Interest¬ 
ingly, the BEC-BCS crossover occurs where color-flavor lock¬ 
ing appears in more elaborate QCD models. The triple point 
(3P) occurs at the baryon mass /x = M, marking the lower 
limit of viability for our model, with the dashed curve extrap¬ 
olation for jjL < M indicating the region where the associated 
inequalities (discussed in Sec. Ill of the main text) cannot be 
satisfied. Note that mo is the constituent quark mass. Sym¬ 
metries of the Eermi surface Z 2 , Z 2 CZ 2 , and Z 4 are indicated 
to connect to the corresponding regions in Eig.[^ /x < M+mo 
corresponds to the region m > me in Eig. ^and /x > M + mo 
to m < TTXc. Directions of vanishing and increasing running 
quark mass m, Z 4 -chiral symmetry restoration, and decon¬ 
fined chiral symmetric (x^) quark-gluon plasma (QGP) are 
shown. 
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the chiral phase boundaries for the Fermi surface 

= ^/{l/m2){mo-^^ + M), (42) 

rPi^) = V(l/m 2 )(mo+M-M), (43) 


where we have set the bare potential = M and taken 
= Sep/i'ip'^ip)- The constraint on conies from the 
lower viable limit for our model: the Fermi surface cannot 
form below the baryon mass and the four inequalities in 
Sec. Ill cannot be satisfied for /i below the triple point. 
The second condition equates g^/3 with the prefactor 
of the density in ej?, a condition required in order to 
force the density dependence in Eqs. ([4^-([4^ to reside 
exclusively in /i. Equations ([4^-([4^ give the boundary 
curves in the T — ji plane plotted in Eig.|^ where we have 
indicated the crossover region by a dashed curve towards 
the lower right corner of the plot. Note that the quark 
chemical potential is directly related to the quark density 
and baryon chemical potential hy g = fiB/3 ^ Uq. 

The resulting T — ji phase diagram at intermediate to 
high densities shown in Eig. is remarkably consistent 
with results obtained using more traditional approaches 
(see, for example, Eig. 1 and discussions in Ref. my 
A key difference in our result is the appearance of a 
quantum BKT transition between the meson dominated 
Hqq) 7 ^ 0 ) the color superconducting {{qq) 7 ^ 0 ) 
regimes. At zero temperature this is driven by decreasing 
the quark chemical potential below a critical value equal 
to the baryon mass M plus the constituent quark mass 
jji = M ^ rriQ. Einally, note that the triple point for the 
deconfined, hadronic, and superconducting phases occurs 
at /i = M. 


VII. CONCLUSION 

We have found that a quantum Berezinskii-Kosterlitz- 
Thouless phase transition separates the superconduct¬ 
ing phase from the ordinary hadronic phase, within the 
Gross-Neveu model of planar QCD modified by an ad¬ 
ditional repulsive vector interaction. In particular, the 
transition is driven by reducing the quark chemical po¬ 
tential /i through a quantum critical point near the 
baryon mass M. At larger values of the chemical po¬ 
tential, the superconducting phase undergoes a BEC- 
BCS crossover in the region near the color-flavor lock¬ 
ing transition that appears in more sophisticated models. 
The topological degrees of freedom associated with the 
QBKT unbinding mechanism are composite quarks near 
the Eermi surface, which we identify as excitations in the 
baryon skin that extends down from the Eermi surface. 
Significantly, our model reproduces the QCD phase di¬ 


agram around the hadron-quark mixed phase, valid for 
/i > M, obtained purely by analysis of Eermi surface 
modes, which we have characterized in extensive detail. 

Several parallels are worth noting with regards to holo¬ 
graphic BKT transitions, beyond the similarities in driv¬ 
ing mechanisms expounded upon particularly in Sec. |V| 
of the present work. Eirst, we recall that at T = 0 the 
chiral field 7 associated with the third Pauli matrix (73 
in Eq. 0 is inactive for g ^ /ic? which case the 
active modes are consistent with the number of spinor 
degrees of freedom in the 3D ^ 2 D projection SU( 2 ) 
^ U(l) C Spin( 2 ). When /i = /ic, however, the van¬ 
ishing mass opens up a large time-like (in the 2D 
theory) length scale ^ I/P 7 , recovering SU( 2 ) symme¬ 
try through the now active generator ( 73 . This picture 
at criticality is highly suggestive of the holographic-type 
BKT phase transitions [19], where the emerging third 
dimension in our problem is reminiscent of the radial di¬ 
rection in anti-de Sitter space that encodes the energy 
scale of the boundary field theory. 

Another interesting point relates to Efimov states and 
their role in holographic BKT transitions. Such three- 
body bound states are expected to play a fundamental 
part in the holographic mechanism by way of an infi¬ 
nite tower of excitations near the two-body dissociation 
threshold [23]. It is curious that the Cooper pairs in 
our model are realized as Eermi surface bound states of 
baryon like modes, states resembling bosonic trimers that 
dissociate in the absence of one constituent boson. This 
similarity appears even more striking by the fact that our 
baryon like modes are in fact composite quarks for which 
our model predicts an infinite tower of bound states re¬ 
lating to the number of vortices attached to an elemen¬ 
tary quark. More generally, in light of the conclusions 
drawn from our results in Sec. [TV| wherein compositeness 
emerges as an inevitable consequence of Dirac particles at 
high densities, we conjecture that a deeper relationship 
must exist between composite fermions in conventional 
condensed matter quantum Hall systems and the math¬ 
ematical structure of Dirac spinors. 

Confirming connections between holographic and four- 
fermion theories demands a rigorous construction of the 
dual gravitational description for our model. Erom a 
practical standpoint, though, this is somewhat academic 
as the primary appeal of holographic techniques is in cir¬ 
cumventing the challenge of solving the boundary field 
theory directly. Indeed, the appeal of the present ap¬ 
proach resides in the notion that a fortuitous choice of 
the correct degrees of freedom and associated interac¬ 
tions in the boundary theory provides the most insight 
into inherently non-perturbative problems. 
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